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This paper analyzes theoretically the signal propagation in spin transport by modulating the 
current passing through magnetic multilayers. Using a macroscopic description of spin transport 
based on the dynamical Boltzmann equation, we show that time-dependent spin transport possesses 
a wave-like character that leads to modifications of pure spin-diffusion dynamics. In particular, the 
wave-like characteristics allow one to extract a finite spin signal-propagation velocity. 
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I. INTRODUCTION 

Time-dependent spin transport in magnetic multilay- 
ers with current perpendicular to the plane (CPP) is 
studied because of its significance in physics and promis- 
ing applications in spintronics devices*^ Most theoreti- 
cal investigations are based on a diffusion equation for the 
spin accumulation or magnetization.—^^ These theories 
show that if one drives a spin-polarized current through 
an interface from a magnetic to a non-magnetic metallic 
layer, the spin propagates by "diffusing" into the non- 
magnetic layers. If one considers time-dependent spin 
transport, such as spin transfer torque switching al- 
ternating current (AC)», or magnetization switching, 6 
where a time-dependent signal is encoded in the spin ori- 
entation, one faces a difficulty of the diffusion equation 
in that no propagation velocity for the spin signal in the 
nonmagnetic layer can be defined. Or, stated differently, 
the diffusion equation yields an infinite propagation ve- 
locity for the spin signal in the metal, because the signal 
will appear everywhere as soon as the source is switched 
on£ In this paper, we show how a physical propagation 
velocity for spin signals in the CPP configuration can be 
determined by deriving and analyzing macroscopic dy- 
namical equations for spin transport. 

It is an interesting connection that a problem anal- 
ogous to that of an infinite signal propagation velocity 
in the spin diffusion equation exists for the heat diffu- 
sion equation, which yields an infinite heat-conduction 
velocity. This difficulty was resolved by recognizing that 
the theoretical description of heat transport needs to be 
generalized by substituting the Maxwell-Cattaneo equa- 
tion 7 -^ for Fourier's law. In this way, one obtains the 
physical picture that heat conduction is characterized 
by a wave-diffusion duality. Formally, the heat diffusion 
equation needs to be replaced by an equation that is es- 
sentially a telegraph equation3iiS As we show in this pa- 
per, a similar modification of the spin-diffusion equation 
is necessary in the case of spin transport. 

We base our derivation of the macroscopic equations 
for spin transport through multilayers on the theory de- 
veloped for steady-state spin transport across magnetic 
multilayers by Valet and Fert^i Instead of using the 
time-independent Boltzmann equation as in Ref. I 111 , we 



treat time-dependent spin transport starting from the 
dynamical Boltzmann equation, which allows us to de- 
rive macroscopic equations and to generalize the spin- 
diffusion equation. 

This paper is organized as follows. The macroscopic 
dynamical equations are derived in the Sec. [IT] of our pa- 
per. Since the central equations (fT5)) and (fll?|) , can also 
be cast in a form that resembles telegraph equations, we 
discuss qualitative aspects of dynamical spin transport 
in Sec. IIIII using these telegraph equations. In Sec. IIV1 
we analyze two concrete examples of time-dependent spin 
transport numerically, and the main conclusions are sum- 
marized in Sec. [V] 

II. TIME-DEPENDENT EQUATION SYSTEM 

In this section, the model of Valet and Fertii for spin- 
dependent transport of conduction electrons through 
metallic multilayers will be extended to take into account 
the time-dependence of spin transport. The electron dis- 
tribution function f s (z, v, t) satisfies the linearized Boltz- 
mann equation 

at oz oe 

= fd 3 v'S[e(v') - s(v)]P s [z, e(v)} [f. (z, v', t) - f.(z, v, t)] 

+ fd 3 v'S[s(v')-e(v)]P s{ [z,e(v)] [/-.(z,v',t)-/ a (*,v,t)], 

(1) 

where — e and e(v) = mv 2 /2 denote, respectively, the 
charge and kinetic energy of the electrons, and E(z, t) = 
—dV(z,t)/dz is the local electric field^ P s (z,e) and 
P s t(z,s) are the spin conserving and spin-flip transition 
probabilities, respectively. Following Ref. [ll|, we assume 
/ s (z,v,t) to be the sum of the Fermi-Dirac distribution 
f°(v) and small perturbations: 

f s (z, v, t) =/°( u )+JL{ [fj,°-^(z, t)] +g s (z, v, t)} , (2) 

where fx — mvp/2 and fx s (z,t) are the equilibrium and 
nonequilibrium chemical potentials, respectively. Due to 
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the cylindrical symmetry of the system around the z axis, 
g s (z, v, t) can be expanded in Legendre polynomials of 
cos 9, where 9 is the angle between v and the z axis, as 



g s (z,v,t) = Y,9i n) (z,t)P n (cos9). (3) 



Here, the zero-order (isotropic) term is absent be- 
cause (df°/de)g s (z,v,i) was defined by Eq. © as the 
anisotropic part of the electron distribution perturbation. 
Using Eq. ([3]), we obtain 



dg s (z,v,t) dg s (z,v,t) 

7u +Vz a 

at dz 



1 1 



dfi s (z,t) dn a (z,t) fi s (z,t) - M-sCM) 

- Wz 



(4) 



dt 



dz 



where p s (z, t) — fi s (z, t) — eV(z, t) is the electrochemical 
potential for electrons with spin s. The derivation of this 
equation is detailed in Appendix IA 21 Note that v in 
Eq. has been restricted to the Fermi velocity vp, i.e., 
|v| = v-p and v z = v-p cos 9. 

With the relaxation times r s and r s f (see Eqs. (|A9|) and 
(|A10[) in Appendix IA 2[) , the local electron mean free path 
X s , diffusion constant D s , and spin diffusion length l s can 
be defined, respectively, as \ s — Vpr' s , D s = ufA s /3, and 
l s = {DgTsf) 1 / 2 , where the momentum relaxation time t' s 
is defined by 



l/r' s = l/r s + l/r si 



(5) 



The appropriate "average" spin-diffusion length l s f can 
be defined as (l// sf ) 2 = {l/l+f + (l/l-) 2 - Throughout 
the paper, subscripts + and — stand for the absolute 
spin directions "up" and "down" , respectively, whereas 
subscripts j and J, stand for the majority and minority 
spin directions, respectively. 

Using the method of Appendix B in Ref . [llj, we express 
the time-dependent current density for spin s as 



J s (z,t) 



e 
V 



^2f s (z,v,t)v z = K 9 P(z,t), (6) 



where k — o~ s /(e\ s ). Note that k is independent of s 
and of the material in the Valet-Fert model. The con- 
ductivity o s can be written as a s = e 2 n s T' s /m, where 
n s = 4ir(mvp /h) 3 /3 is the number of electrons with spin 
s. It is easy to see that a s satisfies Einstein's relation 
o~ s = e 2 N s D s , where 



1 

4^2 



Ns = __( 2m /n 2 ) 3/2 .-" 



/'■'■ (7) 

is the density of states for spin s at the Fermi level /jP, 
and N + = . 

Substituting Eqs. © and © into Eq. (gj), we obtain 

e dJ s {z,t) 1 dfi s (z,t) _ p s (z,t) — p,- a (z,t) 



dz 



ot 



j (8) 



g s d/j, s (z,t) , dJ s (z,t) 
J s {z,t) = — — t s — . (9) 



In steady state, J s (z,i) and fl s (z,t) become time- 
independent and then Eqs. ([8]) and ([9]) reduce to the 
Eqs. (10) and (11) in Ref. fill, respectively. 

Equations © and © will be transformed to more 
directly usable forms next. Without loss of generality, 
the magnetization of the ferromagnet is set to be "up" . 
Then, the majority (minority) spins, which are antipar- 
allel (parallel) to the local magnetization (electron mag- 
netic moment is fi = — (e/m)s) and denoted by subscript 
| (I), point to the absolute spin direction "down" ("up") 
denoted by subscript — (+). In terms of J m (z,t) — 
J + (z,t) — J_(z,t) and (j, m (z,t) = /j,+ (z,t) — /i_(z,i), 
Eqs. ([8]) and ([9]) can be written as 



dJ m (z,t) ^ dfM m (z,t) _ ^ fJ, m (z,t) 



dz 



dt 



(10) 



Jm{z, t) = eN s D t — pJ{z,t), (11) 



dz 



dt 



where 



T x = r sf /2 



(12) 



can be regarded as the spin relaxation timeJ^ The "av- 
erage" diffusion constant D is defined as D = c 2 r, with 
the wavefront velocity c defined by 

c 2 = v 2 /3. (13) 

The "average" momentum relaxation time r is 



l/r=(l/ri + l/ri)/2, 



(14) 



and we have the identity l s i — c\J tT\ . In Eq. (jlip . 
(3 = [t'_ — t' + )/{t'_ + t' + ) equals f3 and for the ferro- 
magnetic and nonmagnetic layers, respectively. The bulk 
spin asymmetry coefficient (3 in the ferromagnetic layer 
is defined by p T(i) = l/cr T(i) = 2p* F [1 - (+)/3], where p F 
is the total resistivity of the ferromagnetic layer. In the 
nonmagnetic layer, we have Pm) = where pjj is the 
total resistivity of the nonmagnetic layer. 

In Eq. pT|) . J(z,t) stands for the total current density 
J(z,t) = J + {z,t) + J-(z, t). By introducing p(z,t) — 
t) + p,— (z,t)] /2, we can also derive equations for 
the charge dynamics 



dJ(z,t) 



dz 



dt 



d 



J{z, t) = 2eN s D— [p{z, t) - eV(z, t)} 
dz 

_ r a7^t)_^ 

dt 



(15) 



(16) 



To describe spin accumulation by spin density instead 
of the chemical potential, it is necessary to transform 
Eqs. dTUJ) an d (fi~Tj) using the following identity (Eq. (|A15j) 
in Appendix IA 3[) 



dz 



dt 



n m (z,t) = ~eN s p m (z,t), 



(17) 
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where n m (z, t) = n+(2, t) — n_(z,t) is the spin density 
and n s (z, t) the nonequilibrium charge density for spin s. 
Using Eq. (fTT|) . we can rewrite Eqs. fTO]) and (fTTjl as 



dJ m (z,t) dn m (z,t) n m (z,t) 



dz 



dt 



(18) 



t / ,n f^ dnmjz^) dJ m (z,t) ~ 
J m (z,t) = -D — — t — pJ{z,t). (19) 
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To proceed further, we need the following identity 
(Eq. (|A16|) in Appendix IA3]) 



n(z, t) - 2n° s = -eN s [2/x(z, t) - 2/x°] , (20) 

where n(z,t) = n-\-(z,t)+n—(z,t) is total nonequilibrium 
charge density and n? s the equilibrium charge density for 
spin s. Using Eq. (|2H)) . we can rewrite Eqs. (fl~5)) and fH 

as 



dJ(z,t) dn(z,t) 



dz 



= o, 



(21) 
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FIG. 1: Variation of Zd/A s , z and A/A s , z with for three 
different values of ry. The short-dashed, solid, and long-dashed 
curves correspond to r\ = 0.1, 0.02, and 0.01, respectively. 



J(z,t) 



dz 
dJ(z,t) 
dt 



dz 



(22) 



/3J m (z,t). 



In general, Eqs. (fT8|) and (|19| should be solved together 
with Eqs. (|2"Tj) . (f2"2"| . and Poisson's equation. However, 
in metals and degenerate semiconductors, the accumula- 
tion of charge occurs on a much smaller length scale and 
varies much faster than that of spini^i^ Thus as an ap- 
proximation, it is assumed that the charge accumulation 
described by n(z, t) can always reach its steady state in- 
stantaneously when spin transport is considered. This 
means that we always set dn(z,t)/dt = 0, which leads 
to dJ(z, t)/dz = according to Eq. (|2~lj) . Therefore, the 
current density J(z, t) in Eqs. (fTTj) and (fT9|) becomes in- 
dependent of z and can be written as J(t) instead. 



III. "TELEGRAPH" EQUATION 

In order to see the physical significance of the dynam- 
ics described by Eqs. (fT5)) and (|19p and to compare it 
with the spin diffusion equation used in Refs. 3.4.5,6, we 
combine Eqs. (fT5)) and HH) to yield the following equa- 
tions 



d 2 n m (z,t) 1 
dt 2 V 



d 2 J m {z,t) 1 

at 2 v 



CM) 



1 <9n 
TV 9* 



J_^9J m (z,t) 



2 d 2 n m (z,t) 
dz 2 



(23) 



2 d 2 J m (z,t) 



dz 2 



dt 



J m (z,t) 
1 dJ{t) J(t) 



dt 



TTi 



(24) 



Because of the formal similarity of each of Eqs. (f2"3"]) and 
(|24| with the telegraph equation, we will refer to them 
as telegraph equations in the following. 

Each of the telegraph equations contains a second- 
order time derivative, which is absent in the spin diffusion 
equation. This term originates from the time derivative 
of the spin current in Eq. (|19p . which is also absent in 
the corresponding equation for the spin current in spin 
diffusion theory; see, for instance, Eq. (8) of Ref. y. This 
additional term shows that it takes a finite time for the 
spin current to adjust to the gradient of the spin accumu- 
lation. 2ils The second-order time and space derivatives 
lead to a wave character of dynamical spin transport in 
addition to its diffusion character described by the first- 
order time and second-order space derivatives. Thus, 
these equations show that time-dependent spin transport 
should be understood using a wave-diffusion duality pic- 
ture. The occurrence of spin accumulation waves enables 
one to determine a well-defined propagation velocity c for 
the signal in time-dependent spin transport. Although 
the spin diffusion equation does not yield spin accumu- 
lation waves and thus a finite signal propagation veloc- 
ity, it can be regarded as an approximation of the wave- 
diffusion duality of the time-dependent spin transport in 
the long-time limit. 

In the following, the telegraph equation of the non- 
magnetic layer will be analyzed in detail. Here, we have 
r = t' s and = 0. Thus, Eqs. ^ and ^ have the 
same structure and we discuss only Eq. (|2"3"]) without loss 
of generality. We seek a damped and dispersive wave 
solution to Eq. (f23|) of the form 



(25) 



n m (z, t) = exp[i(kz — ut)]. 



At this stage, we can set either w = u) T +iu)i or k = k r +iki. 
The complex u> and k will yield damping factors in time 
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and space, respectively. Since we are more interested 
in the damping length (or the dynamical spin diffusion 
length) , we will follow the method of Ref . LL5| and assume 
k = fc r + iki. Substituting Eq. into Eq. (J23J), we get 
the dispersion relation 



iollo + £ = — c k 



(26) 



where a = 1/t' s + 1/Ti and £ = 1/(t' s T\). Separating the 
real and imaginary parts of Eq. (|26|) . we obtain 



1 

2c2 



V(u 2 - 2 + a 2 to 2 ± (lo 2 - 



(27) 



The wavelength, defined as 27r/|fc r |, can be written as 



A 



y/(uj 2 - rf) 2 + {rj + l) 2 Co 2 + (Co 2 - 77) 



where A s . z = ct' s is the z component of the electron mean 
free path. Moreover, we have introduced dimensionless 
quantities 



Co = u)t' s an d V = t' s /Ti. 

The damping length, defined as Id = VI ^i| 
written as 



(29) 
can be 



^(Co 2 - 77)2 + (r, + l) 2 Co 2 - (Co 2 - V ) 



Note that Id can also be regarded as the dynamical spin 
diffusion length. When Co — > or 00, the damping length 
Id will approach Z s f or 21 s {Wt s ^\/(t' s + T\), respectively. 

Figure [1] shows the variation of Zd/A S;Z and A/A SjZ with 
lot' s f° r three different values of 77. Note that the curves 
of \/\ s ,z for different r\ are very close to each other in 
the frequency range shown in the figure. The damping 
length Zd decreases with frequency, which is analogous to 
the skin effect of the electromagnetic wave propagating 
in metal. The intersection of Zd/A S)Z and A/A s z indi- 
cates the critical angular frequency lo c , which separates 
the wave- like region from the diffusion dominated regime, 
because the wave character becomes significant only if 
the damping length exceeds the wavelength. Stated dif- 
ferently, the wave character is significant if the typical 
time scale T s ; g of the time-dependent process is smaller 
than the critical period T c = 2~k/lo c . On the contrary, 
the diffusion character is dominant if r s i g > T c , and the 
spin-diffusion picture becomes a good approximation of 
the wave-diffusion duality in the limit r S i g 3> T c . 

An explicit expression for the critical angular fre- 
quency lo c is obtained by combining A = Id with Eq. (|27[) 



1 r 



7(1 + V) + V7 2 (l + V) 2 +4?? 



~ 7+ (7+ ~)V, 

7 

where 7 = ir — \/(An) 1=3 3.06. 
3.06 + 3.4?y approximately. 



(31) 



Then, we have lo c t' s 




FIG. 2: The variation of v p /c with ujt' s for three different 
values of r\. The short-dashed, solid, and long-dashed curves 
correspond to 77 = 0.1, 0.02, and 0.01, respectively. The thick 
vertical dot-dashed line indicates the critical angular frequen- 
cies uJcTg for the three different values of 77, which are very 
close to each other according to Eq. (|31[). 



The phase velocity, defined as v p = Lo/\k r \, of the spin 
accumulation wave can be written as 



c 



V2 

1 



n - 



\/{Co 2 - rff + + 1) 2 Cj 2 - (lo 2 - i]) 



When Co — > or 00, the phase velocity v p will ap- 
proach 2c/(r; 1 /2 + r; _1 /2) or c, respectively. When 77 = 1, 
the phase velocity becomes equal to c for all frequen- 
cies. Furthermore, the group velocity can be defined as 
v s = dio/dk v and calculated from Eq. (|26[) . 

Figure [2] shows the phase velocity v p as functions of 
lo for 77 = 0.1, 0.02, and 0.01. The phase velocity is ap- 
proximately equal to the wavefront velocity c when the 
wave character is significant (lo > lo c ). In this case, the 
phase velocity provides a good description of the wave- 
like dynamics. On the contrary, when the wave char- 
acter is insignificant (to < lo c ), the wave amplitude is 
damped strongly and the phase velocity is not meaning- 
ful any more. In this region, the propagation velocity is 
the wavefront velocity c, albeit only on the length scale 
of a damping length. 

In the special case where 77 = 1 (t' s = Tx), we have 
\k r \ = w/c and |fci| = l/l s f- This means that the spin 
accumulation wave becomes a non-dispersive but dissipa- 
tive wave with the constant phase (and group) velocity 
c and penetration depth Z s f. However, this case is likely 
not realized because T\ is usually much larger than t' s 
and Valet-Fert theory is justified to be valid only when 
(r^/2T 1 )V2 < L 
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TABLE I: Parameters for Cu and Co used in numerical cal- 
culation. The units of vp, Pn(f)i an d l^f^ are nm/ps, fi-nm 
and nm, respectively, r and Ti are given in ps. 



Material 


Vf 


Pn(f) 


;N(F) 


P 


T 


Ti 


Cu 


1570" 


6 b 


450" 





0.07 e 


3.5 e 


Co 


1570 a 


86 r 


60 d 


0.5 C 


0.005 c 


0.9 e 



"From Ref. |l6t Cu and Co are assumed to have a common Fermi 
velocity in the Valet-Fert model. 
''From Ref. fijl 
c From Ref. [lj 
d From Ref. 03 

^Calculated from Eqs. |fl2t and JT3). 
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FIG. 3: Spin current density J m (z,i) as a function of z. The 
solid curves in (a) and (b) are J m (z, i) at t — 1.75 T a and 
t — 5.75 T a (charge current J(t) = —Jo) with AC drive, re- 
spectively. The dashed curves in (a) and (b) are the spin 
current density J m (z) resulting from the DC current density, 
J = - Jo- 



IV. NUMERICAL RESULTS 

In this section, the general analysis of the telegraph 
equations for spin-transport is augmented by numerical 
solutions for two illustrative examples of signal prop- 
agation using spin polarized currents through a ferro- 
magnet/metal junction: (i) injection of an alternating 
current, and (ii) instantaneous magnetization switching. 
The results are obtained by numerically solving the sys- 
tem of Eqs. (fTU)l and (fTTj) . Our numerical method is 
outlined in Appendix IA 41 These equations are equiva- 
lent to the telegraph equations Eqs. (|23[) and (pM)) . which 
have been discussed in the previous section, but are eas- 
ier to solve. Alternatively, we could solve the equation 




-2000 -1600 -1200 -800 -400 400 
z (nm) 

FIG. 4: Spin accumulation /i m (z,t) for the same parameters 
as in Fig. |3jb). 



system consisting of Eqs. (|T5|) and (fTTj|) . in which the spin 
accumulation is described by the spin density. However, 
it is more convenient to work with the electrochemical 
potential than the spin density when we deal with the 
boundary conditions^ 

We choose a ferromagnet/metal junction consisting of 
Co and Cu as the material system in both of the sce- 
narios. The interface of the junction is placed at z = 
and the Co (Cu) occupies the half-space z > (z < 0). 
The positive direction of the current is parallel to the 
positive direction of the z axis. For simplicity, the inter- 
face resistance of the junction will be neglected. Then, 
the electrochemical potential and the current density are 
continuous across the interface. Consequently, the spin 
transport across two layers can be described by one com- 
mon equation system with different material parameters 
for the two layers. 

The material parameters used in our numerical calcu- 
lation are shown in Table [U All other parameters can 
be obtained from the values in Tab. HI In particular, the 
wavefront velocity is calculated to be c = 910nm/ps from 
Eq. (fT3|) . In the nonmagnetic layer, r] = t' s /T\ = 0.02. 
The wavelength A and damping length are shown as 
the solid curves in Fig. [1] The critical period T c can be 
estimated to be 2t' s ~ 0.14ps from Eq. (f^Tj) . and the 
phase velocity is plotted in Fig. [2] 



A. AC current injection 

The alternating charge current density passing through 
the ferromagnet/metal junction is assumed to be of the 
form J(t) = JosinfW), where Jo = 100 nA/nm 2 . Note 
that the z-dependence of the charge current J(z, t) in 
Eq. (fTTj) is neglected for the investigation of the spin 
transport as pointed out in Sec. [TXJ Two typical fre- 
quencies are studied in the case of the AC drive: v a = 
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FIG. 5: Spin current density J m (z,t) as a function of z. The 
solid curve is J m (z,t) at t = 1.75 Tb (charge current J(t) = 
—Jo) with AC drive. The dashed curve is the spin current 
density J m (z) for the case of a DC current density J = — Jo. 
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FIG. 6: Spin accumulation /i m (z,i) for the same parameters 
as in Fig. [5] 



u) a /(2ir) = 8.33 THz and v h = w fc /(27r) = 0.23 THz, 
which are larger and smaller than the critical frequency 
v c = lo c /(2h) = 7.11 THz of Cu, respectively. The corre- 
sponding periods of the two frequencies are T a = 0.12 ps 
and T b = 4.4ps, which satisfy T a < T c < T b . The nu- 
merical results for the two frequencies are discussed in 
the following. 

High-frequency case (lu — uo a > lu c ): Figure [3] shows 
snapshots of the spin current density J m (z,t) at times 
t = 1.75 T a and t = 5.75 T a . At both times, the charge 
current density J(t) reaches its minimum J(t) = — Jo- 
The wavefront, i.e., the spin signal, can be seen clearly 
in Fig. [3fa), where the time t is so small that the wave- 
front has not propagated beyond the scale of the damping 
length Id. In Fig. Gl^b), the signal has propagated fur- 



ther, and due to the attenuation of the wave, the wave 
front is less clearly visible. Nevertheless, the wavefront 
velocity c can be determined numerically (or experimen- 
tally) by tracking the motion of the wavefront over a 
short time interval after switching on the drive current. 
Since we are using a signal time scale shorter than the 
critical time, we expect from the analysis in Sec. Mil (see 
also Fig. \2§ the phase velocity to be v p ~ c =910nm/ps 
from Eq. (|13p . and a wavelength A = 108 nm. These ex- 
pectations are borne out by the numerical results. The 
dynamical damping length 1^ can also be extracted from 
the numerical data, or from an experiment, by fitting a 
decay time to the envelope of the spin-current signal for 
longer times. Due to inaccuracies of the fitting procedure, 
this quantity is more difficult to determine quantitatively, 
but agrees well with the damping length — 126 nm ex- 
pected from Eq (|30[) . An important qualitative conclu- 
sion can be drawn by comparing the decay of the dynam- 
ical spin signal in Fig. [3Jb) with the spin current density 
J m (z) that results from a DC current density J = — Jo, 
which is also shown. Since our dynamical equations and 
the spin-diffusion equation have the same long-time limit, 
the DC result is identical with steady-state spin diffusion. 
It is apparent that the damping length becomes much 
shorter than the spin diffusion length l s f of the steady- 
state spin transport with DC bias. This is the "skin" 
effect, which is already present in the analytical results 
in Sec. HTT1 

Figure 3] shows the z-dependent spin accumulation /z m 
for the same parameters as in Fig.[3jb). The wavelength, 
damping length, and phase velocity given by Fig. [?] are 
very similar to those in Fig.[3]Jb). Note, however, that the 
amplitude of the dynamical spin accumulation is much 
smaller than the spin accumulation of the steady-state 
spin transport shown by the dashed curve. The rea- 
son is that the AC drive oscillates too fast so that the 
spin accumulation does not have enough time to reach 
its steady-state value. 

Low-frequency case (ui = Wb < oj c ): Before analyz- 
ing the signal propagation velocity, we first show how 
the results change qualitatively compared with the high- 
frequency case. In Figure the spin current J m is plot- 
ted as a function of z driven by an AC current with 
frequency a->b, which is smaller than the critical angular 
frequency oj c . The period Tb of the AC drive is 4.4 ps, 
which is much larger than T a in Fig. [3J The solid curve 
is J m (z,t) at t — 1.75Tb (charge current J = — Jo) with 
AC drive. The dashed curve is again J m (z) driven by 
a DC current density J = —Jo. For this driving fre- 
quency, the wave character is insignificant, because the 
wavelength A = 1856 nm becomes much larger than the 
damping length — 268 nm, and thus the wave ampli- 
tude is damped to zero within just one wavelength. From 
a practical point of view, the wavelength and the phase 
velocity v p — 0.47c lose their meaning in this case. Com- 
parison between the solid and dashed curves shows that 
the damping length for Tb becomes longer than that for 
T a in Fig. O which is a consequence of the "skin" effect. 
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FIG. 7: Spin current density J m (z, t) as a function of z. 
The solid, dot-dashed, and dashed curves are J m (z,t) at 
t = Tb/16, Tb/8, T b /4, respectively, with AC drive. 
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FIG. 8: Spin accumulation fi m (z,t) for the same parameters 
as in Fig. [7] 
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FIG. 9: Spin current density J m (z, t) as a function of z. 
The solid curves in Figs, (a), (b) and (c) are Jm(z,t) at 
t — 0, 0.5, 3.5 ps, respectively. The dashed curves in (b) 
and (c) are J m (z,t) at t = ps plotted again as a reference. 



"wave front" velocity c, which according to our analysis 
of the telegraph equation is still the propagation 
velocity. 

Time-dependent spin transport in the low-frequency 
case can be described approximately by the conventional 
spin diffusion equation. However, it is impossible to es- 
timate the signal propagation velocity from conventional 
spin diffusion theory because there is no wavefront in that 
case and the signal appears in infinity once the charge 
current J(t) is switched on&22- 



B. Magnetization switching 



Figure [5] shows the spin accumulation func- 
tion of z. The parameters used are the same as those in 
Fig. The features of the spin accumulation are again 
reminiscent of the spin current in Fig. [5l Note that the 
spin accumulation has become larger compared with the 
AC drive with period T a in Fig. [H This is reasonable 
because the AC drive oscillates more slowly than that in 
Fig. 0] so that the spin accumulation has more time to 
approach its steady-state value. 

At the time the snapshots in Figs. [5] and [6] are taken, 
no "wave front" of the spin current, or signal, can be 
distinguished. To determine the propagation velocity, we 
show J m (z,t) and li m (z,t) at t = Tb/16, Tb/8, and Tb/4 
in Figs. [7] and [51 respectively. By tracking the motion of 
the wavefront with time, we can estimate the propagation 
velocity of the signal. The result is in agreement with the 



The instantaneous switching of the magnetization in 
the ferromagnet, through which the current passes into 
the nonmagnetic metal, provides perhaps the conceptu- 
ally cleanest picture of a spin-switching process. For 
a numerical study of this process, we consider again a 
ferromagnet /metal junction consisting of Co and Cu. 
We assume that the system is in a steady state in the 
presence of the DC drive with a charge current density 
Jo = 100 nA/nm 2 before the magnetization of the fer- 
romagnetic layer is switched from "up" to "down" at 
t = 0. We model the switching as an idealized instan- 
taneous process, and only consider the evolution of the 
spin current density and spin accumulation afterwards. 
The spin-up electrons become the majority, and the spin- 
down electrons become the minority after the instanta- 
neous switching. The conductivities of the majority and 
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FIG. 10: Spin accumulation /j, m (z,t) for the same parameters 
as in Fig. [9] 



V. SUMMARY 

We studied signal propagation in time-dependent spin 
transport through magnetic multilayers using an exten- 
sion of the Valet-Fert theory to time-dependent phenom- 
ena. We established that time-dependent spin transport 
has a wave character in addition to its diffusive char- 
acter, which enabled us to determine the finite propa- 
gation velocity of signals in spin transport, such as AC 
spin injection and magnetization switching. The propa- 
gation velocity is the wavefront velocity c = vf/\/3- The 
wave character is significant if the signal time scale r s i g 
is smaller than a critical time T c . When the wave char- 
acter is significant (r s i g < T c ), the time-dependent spin 
transport should be modeled by the dynamical equations 
introduced in this paper, or, equivalently, the telegraph 
equations. However, pure diffusive spin transport can be 
regarded as an approximation of the wave-diffusion du- 
ality for slow switching times (r s i g ^> T c ). In this limit, 
the spin diffusion equation can be used to study the time- 
dependence of spin transport approximately, but it incor- 
rectly yields an infinite signal-propagation velocity. 
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minority channels are also exchanged by the switching. 
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that the analysis of the telegraph equation still applies. 



We acknowledge financial support from the state of 
Rhcinland-Pfalz through the MATCOR program, and a 
CPU-time grant from the John von Neumann Institut for 
Computing (NIC) at the Forschungszentrum Julich. 



APPENDIX A: IDENTITIES AND DERIVATIONS 



Figures [9"1 and [TTil show the dynamics of the spin current 
and spin accumulation. Starting from the steady-state 
value shown in parts (a), the magnetization is switched 
instantaneously at t = 0. Figs.OHh) and Fig. [T07 b) show 
snapshots 0.5 ps after the switch, when a pronounced 
kink has developed. This kink indicates the leftmost po- 
sition to which the magnetization-switching signal has 
propagated after 0.5 ps. The kink is noticeable only if 
the time t is so small that it does not propagate be- 
yond the length scale of the spin-diffusion length l s { , over 
which the steady-state signal decays. Thus the signal- 
propagation velocity can be estimated roughly by track- 
ing the motion of the kink with time at the early stage 
of the switching. The result is very close to the wave- 
front velocity c = 910 nm/ps calculated from the ana- 
lytical result, Eq. (fT5|) . Moreover, Fig. [^c) shows that 
the spin current density reaches the steady state with 
"down" magnetization on the time scale of the spin relax- 
ation time T\. Since Ti S> t' s , we can consider t = T\ as 

the long-time limit. This behavior is consistent with the 

i — i 

result calculated from the diffusion equation in Ref. [3J, 
so that again the diffusion character of spin transport 
emerges as an approximation of the wave-diffusion char- 
acter in the long-time limit. 



Useful identities 



Several useful identities will be established by the help 
of Eq. © . Multiplying sin 9 and integrating over 9 from 
to 7r on both sides of Eq. , we have 

/ d9sm9g s (z,v,t) = V^ n) (z,<) / d9 sin 6>P n (cos 9). 
Jo n=1 ' Jo 

(Al) 

The right-hand-side (RHS) of Eq. $K)} can be further 
written as 

RHS = jr g i n) (z,t) f duP (u)P n (u). (A2) 

n=l 

Using the orthogonality relation between Legendre poly- 
nomials 

,i 2 

duP n ,(u)P n (u) = - — — <5 nn /, (A3) 
In + 1 

where 8 n ,ri is the usual Kronecker symbol, we obtain 
from Eqs. (jAlj) and (lA2j) 



(A4) 



d9sm9g s (z,v, t) = 0. 
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Eq. (|A4[) further leads to 

E, . Um 3 f ,^ , , 

g s (z,\,t) = —3- / d vg s (z,v,t) 

Vm 3 r** r* 



Then we can write the RHS (the collision terms) of 
Eq. © as 



ft 3 



dip d9sin0 dvv 2 g s (z, v, <) = 0. 
Jo Jo Jo 



(A5) 



2. Derivation of Eq. (g]) 



Following Ref. [TT|, we substitute Eq. (0) into Eq. flU 
and use the following identity 



9/° _ J_df_ _ -5(v - v F ) 
de mv dv mv-p 



(A6) 



dt 



collision 



9s(z,v,t)-- / dd' shx9'g s (z,v',t)\ , 



■-^-P.[z,e(v)] 

de m 

cte m 

df°(v) „ r ... 47TW r . . 

cte m 



ff s (z,v,i)-- / d0' sin flVsO^v', *)],,,. 



Using Eq. (|A4[) , we can write Eq. (|A7|) in the form 
df s (z,v,t) 



dt 



collision 



-TT 1 ^ t z ' £ u 9s (z, v, t) - -^r^Psf [z, e(v)} g s (z, v, t) 

de m os m 

df°(v) „ , ... A-kv , , , 



(A7) 



(A8) 



By introducing the relaxation times 

W) = PsMv)] ^ (A9) 

-4t = Psf[z,e(v)} — , (A10) 

T s f(«) TO 

where the z-dependence of the relaxation times is ne- 
glected within the same layer, we can further write 
Eq. CM} as 



collision 

df°(v) /ll, 

'- + — )g,(z,v,t) (All) 



df a {z,v,t) 



dt 



+ 



de \t s r sf 
df°(v) fi s (z,t)-fi_ s (z,t) 
de T sf 



3. Derivation of Eqs. (flT|) and ((20)) 

Multiplying both sides of Eq. ([2]) by —e/V, summing 
over v, and using Eq. ()A5|) . we obtain 



,(z,t)-n° s = -eN a \pL a (z,t) - ft ], 



where 



n s (z,t) 



e 
V 
e 
V 



^2f s (z,V,t), 

V 

E/» = - 



(A12) 

(A13) 
(A14) 



are the nonequilibrium and equilibrium charge density 
for spin s, respectively. In turn, Eq. (|Al2j) yields 



n m (z,t) = -eN s fi m (z,t), (A15) 
n(z,t)-2n° = -eN s [2fi(z,t)-2fjP], (A16) 



Taking into account the left-hand-side of Eq. ([T]) and in- 
tegrating over v, we can finally derive Eq. ([4]). Note that 

t s (v) and r s f(v) are restricted to the Fermi velocity vp where n m (z,t) = n + (z,t) — n_(z,t) is the spin density 

after the integration over v and then they are simply and n(z, t) — n + {z, t) + n_(z, t) the total nonequilibrium 



written as r s and r s f. 



charge density. 
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4. Numerical solution of Eqs. (|10[) and 



For the numerical solution of Eqs. (fT0|) and (fTTj) we 
use the method of characteristics and Hartree's compu- 
tational form. Following Ref. [2l|, the space z and time 
t are discretized into grids with equal intervals Az and 
At, respectively. The discretized forms of J m (z,t) and 
/i m (z, t) are J™ i and /i™ 4 at ith space point and nth time 
point, respectively. Then, J™^ 1 and (J^{ at (n + l)th 
time point can be calculated by the iteration relations 



-V 



1 



eiV s 



At 

At 
27 



i) 



(A17) 



J, 



2 + — )J" +1 



1 - ^) (^-i + - v^ (J " + jn+1) ' 

(A18) 

for all space points except the two boundary points, 
which should be determined by boundary conditions. 
Here, J™ is the total current density at nth time point. 
Moreover, Az and At are chosen to satisfy the relation 
Az = cAt. Eqs. (|A17P and (|A18|) can be iterated nu- 
merically to yield the results presented in Sec. IIVI In 
the numerical solution, we used the following initial and 
boundary conditions for the AC spin injection and mag- 
netization switching. 

AC spin injection: The initial conditions are /i m (z, t — 
0) = and J m (z,t = 0) = 0. The boundary condition 
for /i m (z,t) is fim(z — ±oo, t) — 0. From Eq. (fTTj) . the 
boundary condition 



Jm{z — ±00, t) 



(3Jo 



1+L0 2 T 



, ,uit cosfwt) — sin(wt) 



(A19) 



-exp(-t/r)]. 



for J m (z,t) can be derived. 

Magnetization switching: The initial conditions for 
J m (z,t) and li m (z, t) are the steady-state solutions to 
Eqs. dTOJ) and (fl"T|) 

= 0) = C exp(-z/; F f ), (A20) 
J F (z,t = 0) = --^^exp(-z/Z F )-/3J ,(A21) 



l£(z,t = 0) = 0>exp(z/Z s N f ) 

jN( 2 ,t = 0) 



N/ . n\ ^0 / /,N\ 



2epU 



(A22) 
(A23) 



N l sf 



where C = -2e^J (p^ s F f ^e f )/(PF's F f + PnO- Here > 
/Lt F and J F apply to the ferromagnetic layer occupying 
z > 0, whereas /u^ and J^J refer to the nonmagnetic layer 
[z < 0). In deriving the initial conditions above, we have 
used the identity 1/(2/?^^) = ctn(f)/2 = e 2 N s Dn(f), 
where <7n(f) is the total conductivity of the nonmag- 
netic (ferromagnetic) layer. The boundary condition for 
fi m (z,t) is n m (z = ±oo, t) = 0. Then, the boundary con- 
dition for J m (z,t) can again be derived from Eq. (jlip . 
This yields 



J m (z = ±oo, t) = PJ [1 - 2exp(-t/r)], 



(A24) 



where /? is the asymmetry parameter before the magne- 
tization switching. Note that (3 becomes —(3 when the 
magnetization is switched (t > 0). 
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